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o 
o 

^ ' For the single-molecule magnet Mni2 , Cornia et al. recently proposed that solvent molecules may 

^ ' cause the quantum tunneling that requires a lower symmetry than S4. However, magnetic quantum 

tunneling and electron paramagnetic resonance experiments suggested that the proposed theory may 
not correspond to the measurements. In this regard, we consider positional disorder induced by the 
' solvent molecules and orientational disorder by the methyl groups of a Mni2 molecule. We calculate, 

within density-functional theory, the second-order transverse magnetic anisotropy parameter E and 
an easy-axis tilting angle 6 induced by the positional disorder and the E value by the orientational 
disorder. We also calculate the local magnetic anisotropy and the local easy axis for each inequivalent 
O , Mn site in different environments to investigate their effects on the global E value. We find that 

. the hydrogen bonding between a Mni2 molecule and the solvent molecule is crucial to obtain a 

substantial E value and that the E value increases upon geometry relaxations. Our calculations on 
relaxed geometries show that the largest calculated E value is 0.016 K and that the largest tilting 
angle is 0.5°. Our largest E value is comparable to experimental results but larger than Cornia et 
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PACS numbers; 75.50.Xx, 75.45.-l-j, 75.30.Gw, 71.15.Mb 



I. INTRODUCTION 



o , 

O , Single-molecule magnets (SMMs) are three-dimensional arrays of identical molecules. Each molecule has a large 
effective spin with a large magnetic anisotropy barrier, and responds to an external magnetic field as a nanoscale, single- 
domain magnetic particle. As such SMMs are ideal for high-density magnetic storage devices. The most extensively 
> . studied SMM is [Mni20i2(CH3COO)i6(H20)4] ■2(CH3COOH)-4(H20) (hereafter Mni2)i A single molecule of Mni2 
^-H ' (Fig.ni has an effective ground-state spin of 5 = 10 with a magnetic easy axis along the crystal c axis (i.e., z axis), and 
has 84 symmetry. Low-temperature magnetization hysteresis loop measurements on Mni2 revealed steps at quantized 
magnetic fields, in which macroscopic quantum tunneling (MQT) occurred between spin-up states and spin-down 
I , states through the magnetic anisotropy barrier. Many competing models have been proposed to understand various 
' features of the MQT: fourth-order transverse anisotropy,^ thermally assisted quantum tunneling,-'^ the Landau-Zener 
effect,^ and dipolar interactions with dynamic hyperfine fieldsii However, the observation of quantized steps that 
"Y^ - are forbidden by the symmetry of a single Mni2 molecule remains a long-standing puzzle. It is imperative to fully 
' understand the MQT because higher tunneling rates imply the loss of magnetically stored information, while tunneling 
^ ,, is useful in quantum superposition of states crucial to quantum computing?.. 
I ' A general effective single-spin Hamiltonian for SMMs is, to fourth order, 
T3 ' 

n = DSl + E{Sl-Sl) + CiSl + C2{SX + St)+n' + giiBB-S, (1) 

' where D is the uniaxial anisotropy parameter, E is the second-order transverse anisotropy parameter, Sz is a spin 
^ : operator projected onto the easy axis, Ci and C2 are the fourfold longitudinal and transverse anisotropy parameters, 
Ti! contains all possible fourth-order terms that break S4 symmetry, and the last term is the Zeeman interaction. 
The experimental parameter values are D = -0.55 K, Ci = -0.00117 K, and C2 = ±2.88x10^5 KiS The MQT is 
d ' ascribed to the transverse terms in Eq. (Q. According to the S4 symmetry of a Mni2 molecule, E = and TC' = 
so that the MQT is caused by C2 alone and is allowed between states whose eigenvalues of Sz differ by integer 
multiples of 4. However MQT is observed between states that violate this rule. To understand these unexplained 
tunneling steps, it is necessary to include a lower symmetry than S4. Transverse fields caused by internal fields 
(exchange, dipolar, and hyperfine fields) or field misalignment could induce these steps, but resulting tunneling rates 
will be too small compared to measured valuesiS. Thus, a nonzero E value is needed to expedite tunneling assuming 
that H' is negligible. So far, two different theories have been proposed in this regard. Chudnovksy and Garanin^i 
proposed that possible dislocations in single crystals may provide a broad continuous distribution of E and showed 
that a broad continuous range of E values is required to explain experiments. Detailed x-ray diffraction studies by 
Cornia et ali^ suggested that disorder caused by acetic-acid (referred to as HAc) solvent molecules may break the 
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S4 symmetry and induce substantial locally varying E values. As discussed in Sec. IIA, such disorder is expected on 
general grounds. Recently, two completely different experimentsi^iii have been performed to elucidate the sources of 
a low symmetry in Mni2. Magnetic quantum tunneling measurementsi^ showed that tunneling splittings (between 
states whose eigenvalues of Sz do not differ by integer multiples of 4) have much narrower distributions than what 
the dislocation-induced theory suggested. They also showed that the largest E value for low symmetry molecules is 
0.01 K, which is twice larger than what Cornia et al}'^ calculated. Electron paramagnetic resonance measurements 
showed that an upper bound on the E values is 0.014 K/^ and that some molecules may have an appreciable easy-axis 
tilt away from the crystal c axis.^^ 

The two different experiments revealed consistent results that may not be fully explained by either of the two 
proposed theories. Thus, there is a need for a new theory or a more refined theory that is consistent with the 
experimental results. In this regard, we investigate the sources of a low symmetry in Mni2 within density-functional 
theory (DFT). Specifically, we consider both positional disorder due to the HAc solvent molecules and orientational 
disorder caused by methyl-group rotations in Mni2 molecules. We find that hydrogen bonding between the Mni2 
molecules and the solvent molecules is crucial in the second-order transverse anisotropy. We calculate the magnetic 
anisotropy parameters and the easy-axis tilting angles from the crystal c axis for representative configurations induced 
by the positional disorder that will be described later. We find that the E value increases upon geometry relaxation, 
and that there exist relationships, due to symmetry, between the E values for the representative configurations. We 
also calculate torsional barriers and the magnetic anisotropy parameters for configurations induced by the orientational 
disorder. Our calculated E values are comparable to the experimental data, while the easy-axis tilting angles remain 
small. In Sec. II, we describe our method and model for the positional disorder and for the orientational disorder. In 
Sec. Ill, we show our results and discuss their consequences. In Sec. IV, we present our conclusion. 

II. METHOD AND MODEL 

In our DFT calculationsi^, we use spin-polarized all-electron Gaussian-orbital-based Naval Research Labora- 
tory Molecular Orbital Library (NRLMOL)iIi within the Perdew-Burke-Ernzerhof (PBE) generalized-gradient ap- 
proximation (GGA) for the exchange-correlation potential. ■^^ We consider the following simplified form of Mni2: 
[Mni20i2(CH3COO)8(HCOO)8(H20)4] • 2(CH3COOH) where 8 formates (HCOO) substituted for 8 acetates 
(CH3COO, referred to as Ac) bridging Mn^+ ions in the inner cubane and Mn'^^ ions in the outer crown (Fig. Q). 
For simplicity, we do not include the water molecules of crystallization. Instead of a periodic structure, we examine 
an isolated Mni2 molecule surrounded by four HAc solvent molecules. The zero-field magnetic anisotropy barrier for 
the S* = 10 ground-state manifold does not change much with our simplification nor with attachment of the solvent 
molecules. The total magnetic moment for the ground state was confirmed to be 20/j,b, which is in good agreement 
with experiment. It is common that the x-ray deduced experimental geometry has underestimated C-H and 0-H bond 
lengths in comparison to standard hydrogen bond lengths, which results in large self-consistent forces on hydrogen 
atoms as large as an order of 1 hartree/bohrjiSi Therefore, we use as an initial geometry the experimental geometry 
with hydrogen positions corrected such that all C-H (0-H) bond lengths become 1.1 A (0.96 A). Then the initial 
geometry is self-consistently relaxed using NRLMOL until forces exerted on all atoms are small. 

A. Positional disorder 

A single Mni2 molecule has four HAc solvent molecules, each of which is shared by two neighboring Mni2 molecules. 
From the perspective of the HAc molecule, there are two equivalent ways to place one HAc molecule between two 
Mni2 molecules. That is, the methyl group (CH3) in HAc can point toward a Mni2 molecule or away from it (Fig.[2Jl. 
In the latter case, there is hydrogen bonding between the oxygen atom 0(2) in the Mni2 molecule and the oxygen 
atom 0(1) in the solvent molecule (Fig.^. We refer to this configuration as head and the other as tail (Fig.[2Jl. There 
are a total number of 16 ways to place four HAc molecules around a Mni2 molecule. Only six of them are inequivalent 
by symmetry as illustrated in Fig. |21 where n denotes the number of the head configurations. The populations of the 
six configurations are shown in Tabled The configurations n — 1, n — 2 trans, n — 2 cis, and n = 3 break the S4 
symmetry, while n — Q and n — 4 preserve the S4 symmetry. The three configurations, n — 1, n — 2 cis, and n — 3, 
can cause easy-axis tilts from the crystal c axis which coincides with the z axis, the easy axis without considering the 
positional disorder. 

To examine if the head configuration affects the magnetic anisotropy in the same way as the tail, we consider the 
following three configurations: (i) two tails only in a twofold symmetric way (trans), (ii) three heads and one tail, 
and (iii) three heads only. We calculate the D and E values for an initial geometry of each configuration. We find 
that configuration (i) provides D = —0.56 K and E = 0.00053 K and that configurations (ii) and (iii) provide the 
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same parameter values oi D = —0.58 K and E = 0.0035 K. These results imply that the tail configuration does not 
substantially contribute to the second-order transverse anisotropy. It is thus sufficient to carry out DFT calculations 
on head-only configurations, which reduces significantly the computation time. Hereafter, unless specified, when we 
refer to the n = 1 configuration, it indicates only one head HAc added to a Mni2 molecule without three tail HAc 
molecules (see the leftmost figure in Fig. O. The same rule is applied to the rest of the configurations. Wc also find 
that both the transverse magnetic anisotropy and the easy-axis tilting angle increase by a factor of 3 upon geometry 
relaxation. 

Especially since we relax an isolated Mni2 molecule with solvent instead of a periodic structure, it is worthwhile 
to check if DFT provides a correct hydrogen bond between two isolated atoms using a simple example system. For 
this purpose, we carry out DFT calculations for a water dimer. Using its fully relaxed geometry, we obtain the 
distance between the two oxygen atoms, the 0-H-O angle, and its binding energy. They are in excellent agreement 
with experimental datc>22i2iiS^ as shown in Table ^1 This result confirms that DFT captures the hydrogen bonding 
between Mni2 molecules and solvent molecules, and that relaxation of the initial geometries should provide physically 
refined geometries. 

With a well relaxed geometry a second-order single-spin Hamiltonian is calculated considering spin-orbit coupling 
only.^'^ The Hamiltonian has a form of ^^^^ ^ ^ ^^vS^S^. Here is obtained from the calculated occupied and 
unoccupied states of the well relaxed geometry and contains the information on the second-order magnetic anisotropy 
parameters and the principal axes of the magnetic anisotropy. By exploiting symmetry arguments appropriate to the 
Mni2 and the solvent molecules, we can explore all of the six different configurations (Fig.|21) using only the following 
two configurations. Once the single-spin Hamiltonians for the n — Q and n = 1 configurations are calculated, those 
for the remaining four configurations can be obtained by rotations of coordinates. For instance, we can describe the 
Hamiltonian for n = 2 trans (cis), Ti.[n = 2], in terms of those for n — 0, Ti.[n — 0], and n = 1, Ti.[n = 1], using a 
rotation matrix R (see Fig.|3}. 

H[n^O] EE J2 ^ISS^S,, n[n^l]^ E , (2) 

fi.,L^—x,y,z fi,iy—x,y,z 

n[n = 2] = E ll!J[S^S. + RiS^)R{S,)]~ ^ 

l_L.v—x.y,z fi,v—x,y,z 

- E 7^1*^'^^^^^.^. (3) 

f_L^i^—x,y,z 

Txy ^ ^T^y^ ^xz ^ ^1 '^yz ^ ^7 ('^) 

(2c) ^ Q (2c)^ (1)_ (1) (2c)^ (1), (1) 
Ixy ^' ixz Ixz lyz i lyz !xz ' lyz 

(3) ^ (1) (3) ^ _ (1) (3) ^ (1) /g^ 
Ixy Ixy 7 Ixz lyz i lyz Ixz \^ / 

where ^^^*\ 7*^^'^\ and 7*^'^^ stand for the 7 matrices for the n = 2 trans, n = 2 cis, and n ~ 3 configurations, 
respectively. Here we show only the off-diagonal elements of the 7 matrices. Since n = has the S4 symmetry, 
Jxy =Jxz ="Yyz =0- Wc have also verified this numerically. 



B. Orientational disorder 



The methyl group in Ac can rotate about the C-C single bond. During the rotation, the energy of the methyl 
group will change. The difference between the maximum and the minimum energy is called a torsional barrier. For an 
isolated Ac molecule, the torsional barrier is calculated to be about 320 K, which is comparable to the thermal energy 
at room temperature. So it is interesting to examine if the orientational disorder caused by rotations of the methyl 
groups can induce the appreciable second-order transverse anisotropy. We consider rotations of the methyl groups 
in a Mni2 molecule, and there are three inequivalent Ac sites. The methyl groups in the Ac molecules bridging the 
inner Mn ions and the outer Mn ions will have less freedom to rotate and therefore a smaller effect on the transverse 
anisotropy. The methyl groups of type A (Fig. are close to the HAc solvent molecules so they also have less 
freedom to rotate. Thus we rotate the methyl groups of type B (in the Ac molecules that are free from the HAc 
molecules) with the methyl groups of type A fixed. In this study, we do not allow solvent-induced symmetry breaking 
and do not consider geometry relaxation induced by methyl rotations. We consider an initial geometry of the n — A 
configuration with several methyl groups of type B rotated, and calculate the torsional barriers and the magnetic 
anisotropy parameters. 
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III. RESULTS 
A. Positional disorder 

The well relaxed geometry for the n = configuration gives rise to D = —0.54 K and E = 0. The well relaxed 
geometry for the n — 1 configuration provides D = —0.54 K, E — 0.008 K, and the easy-axis tilting angle 6 = 0.4°. 
Using this result and Eqs. (0)-®, we calculate D, E, and 9 for the n — 2 trans, n ~ 2 cis, n — 3, and n = 4 
configurations (TableHJ. The D value does not vary much among the different configurations. (Note that D = —0.55 K 
without solvent molecules. The E value for n = f is the same as that for n = 3 and is one half of that for n = 2 
trans. The E value for n = 2 cis is two orders of magnitude smaller that that for n = 1. The largest E value arises 
from n = 2 trans, while the largest comes from n = 2 cis. The E values for the low-symmetry configurations are 
determined by all the off-diagonal elements of the 7*^^-' matrix, while the easy-axis tilts are due only to the nonzero 
ji]} and 7yi^ values [Eqs. 0-®]- For the positional disorder, the magnitudes of all the off-diagonal elements of the 
7^"'^) matrix are very small compared to those of the diagonal elements. Therefore, there is no appreciable easy-axis 
tilt and there is a relation between the E values for different configurations. The n = 2 cis configuration does not have 
a significant contribution to the transverse anisotropy. If the magnitudes of the off-diagonal elements are significant, 
there is a substantial easy-axis tilt and the E values for the distinctive configurations are not correlated. Additionally, 
the E value for n = 2 cis has the same order of magnitude as those for the other low-symmetry configurations. This 
occurs when extra electrons are donated to a Mni2 molecule since the locations of the extra electrons are intrinsically 
dependent on donor locations. In comparison to Cornia et al.'s resultsji? our calculated E values are larger by a 
factor of 3 and our easy-axis tilting angle is almost the same. In comparison to the experimental data^i^ii^ii^ our 
largest E value is in good agreement but our largest easy-axis tilting angle is smaller than 1°. 

To examine the effect of the positional disorder on the local magnetic anisotropy, we calculate the single-ion 
anisotropy parameters'^ and the local magnetic easy-axis tilts for the three inequivalent Mn sites in two different 
environments. These three sites are the cubane Mn^+ ion (at site c in Fig.^ and two Mn'^+ ions (at sites r and b in 
Fig. 0. The two environments correspond to (i) the head and (ii) tail configurations. For the environments (i) and 
(ii) we use a well relaxed n = 2 trans (two heads and two tails) geometry. Before showing our results, we describe 
our method. The lowest-order energy shift A2 due to spin-orbit coupling Vls 
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where Cjo- i^ja-') denotes the energy of the ith occupied (jth unoccupied) unperturbed state whose wavefunction is 
dV'jo-'))- The first sum runs over spin- up and spin-down states and the second sum over all occupied and 
unoccupied states. A wavefunction of the system is expanded in a sum of contributions associated with each atom, 
\'^i<y)='^a IV'S-) ^^'^ li^ja')— J2a I'^fa') ^ whcrc the sums run over all atoms. Thus, A2 can be written in terms of 
projected terms onto each atom and mixed terms between different atoms as follows: 

where the first sum of the mixed term runs over all cases where at least two of the indices aPjS are different. Based on 
the assumption that the contributions associated with each atom, \ipf„), have small interatomic overlap, which is the 
case for Mni2, the mixed terms are expected to be negligible. This is true considering the relative magnitude of the 
Coulomb potential near an atom compared to the interstitial region. Thus A2 is approximately equal to the tensor 
sum of the projected terms onto all atoms, the first term in Eq. JHJ. Regardless of the types of the environments, a 
cubane Mn'*+ ion has the local easy axis almost in the xy plane (the local E value is of the same order as the local D 
value) and the two crown Mn^+ ions have the local easy axis close to the z axis. Therefore the Mn'^+ ions contribute 
most of the total magnetic anisotropy. Since the Mn^+ ions are far away from the solvent molecules, the local E value 
for the head configuration is the same as that for the tail. For the two inequivalent Mn'^+ ions, the local E value is 
sensitive to the environments, while the local D value and the local easy- axis tilting angles are not sensitive to them. 
Especially for the Mn^+ ions of type b, the local E value for the tail configuration is an order of magnitude smaller 
than that for the head. This large difference in the local E values of the head and tail configurations leads to the 
large global E value. 

We also calculate tunnel splitting A_io,4 between Ms = — 10 and Ms — 4, where Ms is an eigenvalue of Sz (Table I) 
for the six distinctive configurations perturbatively and by exact diagonalization. We use our calculated values for D 
and E and the experimental values^ for Ci and C2. The tunneling between Ms = — 10 and Ms = 4 is allowed due to 
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the nonzero second-order transverse anisotropy but the magnitude of the tunnel spUtting is mostly governed by the 
fourth-order transverse anisotropy. So the tunnel splitting is very dependent on the C2 value but less on the E value. 
As shown in Table I, if we reduce the C2 value by one half from the experimental value while keeping other parameter 
values constant, the tunnel splitting is reduced by an order of magnitude. Thus, the experimental value of A„in a[ioI| 
may not be a good indicator of the magnitude of the E value. 

B. Orientational disorder 

When the four methyl groups of type B are rotated relative to those of type A in the n = 4 configuration, we find 
that there are at least four local energy minima at the rotation angles of —10°, —25°, 0°, and 30° where 0° denotes 
the experimental geometry. Since the four methyl groups are rotated simultaneously, we have to divide by 4 to obtain 
the torsional barriers per methyl group for each rotation angle. As shown in Table HVl the torsional barriers at the 
three nonzero rotation angles are comparable to the thermal energy at room temperature. So we consider the three 
twofold symmetric configurations, which correspond to n = 4 with two methyl groups of type B across from each 
other rotated by —10°, —25°, and 30°. Our calculations show that the E values for all three angles are on the order 
of 10~* K, which is much smaller than those due to the positional disorder f Table HV|) . 

IV. CONCLUSION 

Considering the solvent positional disorder and the orientational disorder of the methyl groups in the Mni2 
molecules, we investigated, within density-functional theory, the sources of the second-order transverse magnetic 
anisotropy in Mni2. We found that only hydrogen-bonded solvent molecules contribute to the symmetry-breaking 
transverse anisotropy, and that geometry relaxation is necessary from the standpoint of quantitative prediction. For 
the positional disorder, the large E value is due to the large difference between the local E values of the two different 
environments (the head and tail configurations). There are no significant easy- axis tilts from the z axis for the low- 
symmetry configurations. Accordingly, the E values for the different configurations are correlated. The orientational 
disorder contributes to the E value much less than the positional disorder. Our largest calculated E value from the 
positional disorder quantitatively agrees with the experimental data. However unresolved issues remain. That is, the 
experimental E values suggest a continuous distribution rather than a discrete distribution that is predicted by the 
positional disorder. 
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TABLE I: Populations, calculated magnetic anisotropy parameters {D and E), easy-axis tilting angles 9, and tunnel splittings 
between Ms = — 10 and Ms = 4, A_io,4, for six distinctive configurations of an isolated Mni2 molecule with four HAc solvent 
molecules. See Fig.^l For A_io,4[A], we use Ci = -0.00117 K and C2=2.88x lO"'^ K.^ For A_io,4[B], we use Ci = -0.00117 K 
and C2 = 1.44xl0~'^ K. 



configuration 


population 


D (K) 


E{K) 


9 (degree) 


A_io,4[A] 


A_io,4[B] 


n = 


1 

16 


-0.54 














n = 4 


1 

Ifi 


-0.56 














n = 1 


4 
16 


-0.54 


0.008 


0.4 


7.0 X lO"'^ 


7.8 X 10"* 


71 = 2 cis 


4 
16 


-0.55 


0.00002 


0.5 


1.8 X 10"^ 


4.4 X 10"" 


n = 2 trans 


2 
16 


-0.55 


0.016 





1.0 X lO"'* 


9.3 X 10"* 


n — 3 


4 
16 


-0.55 


0.008 


0.4 


6.7 X lO"'^ 


8.1 X 10"* 



TABLE IL DFT calculated hydrogen bond length between the two oxygen atoms i?(0-0), the bond angle (0-H-O), and its 
binding energy for a water dimer in comparison to experimental datai— 





R(O-O) (A) 


Z (0-H-O) (degree) 


binding energy (eV) 


DFT 


2.906 


174 


0.223 


Exp. 


2.946 [20] 


174 [21] 


0.232 ± 0.030 [22] 



TABLE IIL Calculated local magnetic anisotropy parameters {D and E), local easy-axis tilting angles S, and local anisotropy 
barriers for the three inequivalent Mn sites in different environments. Here "c", "r", and "b" stand for the three sites shown 
in Fig.0 "head" and "tail" denote the head and tail configurations of the HAc solvent molecules. 



site 


D{K) 


E (K) 


S (degree) 


barrier (K) 


Mn''+(c, head) 


-0.50 


0.18 


101.5 


0.72 


Mn''+(c, tail) 


-0.52 


0.18 


88.4 


0.77 


Mn^+(r, head) 


-2.44 


0.037 


11.8 


9.6 


Mn^+(r, tail) 


-2.47 


0.016 


11.4 


9.8 


Mn^+(b, head) 


-2.54 


0.23 


34.5 


9.2 


Mn^+(b, tail) 


-2.65 


0.036 


36.2 


10.4 



TABLE IV: Calculated torsional barriers per methyl group and magnetic anisotropy parameters due to rotations of two methyl 
groups of type B (Fig. [TJ across from each other in the n = 4 configuration. The barriers are calculated from the closest local 
minimum configurations. 



rotation angle (degree) 


torsional barrier (K) 


D (K) 


E{K) 


-10 


316 


-0.57 


0.00009 


-25 


308 


-0.57 


0.0002 


30 


60 


-0.57 


0.0001 
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FIG. 1: A single Mni2 molecule surrounded by four HAc solvent molecules, [Mni20i2(CH3COO)8(HCOO)8(H20)4]4(CH3COOH) 
The sizes of atoms decrease in the order of Mn, O, C, and H. The four inner cubane Mn^^ ions {S = 3/2) are ferromagnetically 
coupled to each other and the eight outer crown Mn='+ {S = 2) ions are ferromagnetically coupled. The Mn spins in the crown 
are antiparallel to those in the cubane, which results in the total effective ground-state spin of 5* = 10. The easy axis of the 
molecule is normal to the page. Two inequivalent acetates are labeled as "A" and "B" . The hydrogen bond is formed between 
oxygen atoms 0(1) and 0(2). 




FIG. 3: Three head-only configurations used for Eqs. Q-©- 



